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Abstract 


Chemical hydrides have been identified as a potential medium for on-board hydrogen storage, one of the most challenging technical barriers 
to the prospective transition from gasoline to hydrogen-powered vehicles. Systematic study of the feasibility of the sodium borohydride systems, 
and chemical-hydride systems more generally, requires detailed kinetic studies of the reaction for use in reactor modeling and system-level 
experiments. This work reports an experimental study of the kinetics of sodium borohydride hydrolysis with a Ru-on-carbon catalyst and a 
Langmuir-Hinshelwood kinetic model developed based on experimental data. The model assumes that the reaction consists of two important steps: 
the equilibrated adsorption of sodium borohydride on the surface of the catalyst and the reaction of the adsorbed species. The model successfully 
captures both the reaction’s zero-order behavior at low temperatures and the first-order behavior at higher temperatures. Reaction rate constants at 
different temperatures are determined from the experimental data, and the activation energy is found to be 66.9 kJ mol~! from an Arrhenius plot. 
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1. Introduction 


The prospect of peaking world-wide oil production has moti- 
vated the development of alternative energy sources and energy 
carriers. Hydrogen has been identified as a likely replacement 
energy carrier, and renewable energy sources such as solar, wind 
and hydro have been identified as potential drivers for future 
energy systems. Producing hydrogen from renewable sources 
and utilizing hydrogen in high-efficiency fuel cells or engines 
may significantly reduce greenhouse gas emissions and pol- 
lution. Furthermore, a hydrogen energy economy would be 
less dependent on limited fossil fuel supplies. Nevertheless, 
many critical technical challenges remain to be addressed before 
hydrogen-based energy can become widely available and eco- 
nomical [1]. 

On-board hydrogen storage has been identified as one of the 
most challenging technical barriers to the transition from gaso- 
line to hydrogen-powered vehicles. Two promising hydrogen 
storage material categories are metal hydrides and chemical 
hydrides, both of which rely on catalysis to improve kinet- 
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ics. The discovery of the catalytic effect of Ti compounds in 
the hydrogenation and dehydrogenation of sodium alanates by 
Bogdanovic and Schwickardi [2] has aroused increased inter- 
est in metal hydride research. The successful demonstration 
of Millenium Cell’s sodium borohydride (NaBH4) system in 
Daimler Chrysler’s fuel cell minivan represents another poten- 
tial hydrogen storage method [3]. Chemical hydride systems 
always involve reactions between the storage medium and a sol- 
vent, in most cases water. The reaction is generally catalyzed, 
and solid catalyst is preferred to facilitate recycling. To study the 
feasibility of sodium borohydride systems and other chemical- 
hydride systems systematically, detailed kinetic studies of the 
reaction, system-level experimental testing, and reactor model- 
ing are essential. 

Much prior catalysis research has focused on gas-phase 
processes primarily used in the petroleum industry. However, 
solid-catalyzed liquid-phase reactions, an area that has not 
received much attention in the past, are crucial for chemical- 
hydride systems. For sodium borohydride systems, the chemical 
hydride is premixed with water to form aqueous NaBH; solution. 
At room temperature, the reaction between sodium borohydride 
and water is very slow. To further suppress the reaction, a small 
amount of sodium hydroxide (NaOH) is commonly added to 
the NaBH; solution. Hydrogen can be liberated by contacting 
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Nomenclature 

Ca concentration of NaBH; at time ¢ (kmol m~?) 

Cao initial concentration of NaBH; (kmol m~?) 

d diameter of the catalyst (m) 

AEsct activation energy for the catalysis process 
(kJ mol!) 

AHads heat of adsorption (kJ mol`!) 

AH,xn heat of reaction (kJ mol`!) 

ko zero-order rate constant (kmol s7! kg catalyst” !) 

kı first-order rate constant (m? s7! kg7!) 

kL rate constant for Langmuir-Hinshelwood kinetic 
model (kmol s~! kg catalyst~!) 

k” normalized rate constant for Langmuir- 


Hinshelwood kinetic model based on reactor 
volume (kmol m~? s~! kg catalyst~!) 


kn nth-order rate constant 
(m?” kmol ~” ) s7! kg catalyst!) 

Kslope slope of the linear line (kmol m7? s7!) 

K Langmuir adsorption isotherm constant 
(m? kmol~!) 

Ko Langmuir adsorption isotherm constant at 25°C 
(m? kmol!) 

Meat weight of catalyst (mg) 

R universal gas constant (8.31 kJ kmol”! K7!) 


AS? entropy change of the adsorption process 
(kJ kmol~! K7!) 


t time (s) 

T temperature (°C) 

To reference temperature, 25°C 

V volume of the batch reactor (m?) 


Vp, the volume of Hz collected at time ¢ (mL) 
VH, max the total volume of H2 collected at the end of 
experiment (mL) 


Greek symbol 
0A coverage of adsorbed BH4~ on catalyst surface 


the aqueous NaBH; solution with a ruthenium (Ru) catalyst, as 
follows: 


NaBH, + 4H20 S$ NaB(OH), + 4H> + Heat (1) 
Ru(cat.) 
The foregoing reaction is highly exothermic, with a heat of 
reaction of 210kJ mol! [4]. Previous researchers have con- 
ducted preliminary studies [5-7] on its kinetics; nevertheless, 
these studies are predominantly qualitative. Furthermore, most 
previous research employed custom-made catalysts, which may 
not be stable and are difficult to replicate. Conversely, commer- 
cially available Ru catalysts offer good stability and consistency. 
Therefore, the main purpose of this paper is to quantify the kinet- 
ics of the liquid-phase reaction over a solid catalyst — sodium 
borohydride hydrolysis catalyzed by commercially available Ru 
catalyst — both to characterize the order of the reaction, and to 
provide a kinetic model for future reactor development. While 
the pH level [8,9] and other catalytic alloys [9,10] are known to 


affect the hydrolysis process, these factors were not considered 
in the present work because: (a) a pH level of 14 is the most 
likely for practical applications involving stabilized solutions 
[11], and (b) Ru remains the most common catalyst for such 
reactions. 


2. Experimental setup 


Accurate study of kinetics requires conducting experiments at 
constant temperature. The NaBHy hydrolysis reaction is highly 
exothermic; therefore, it is essential to maintain constant solu- 
tion temperature by reducing the heat generation rate during 
experiments and by enhancing external cooling. Two methods 
have been used to minimize the heat generation rate during the 
present experiments: the use of a low concentration of NaBH4, 
e.g. 1%, and use of a relatively small amount of catalyst to limit 
the reaction rate to a measurable range. Because the reaction 
rate increases significantly with temperature, the amount of cat- 
alyst used in the experiments was decreased for higher reaction 
temperatures. Type-T (copper—constantan) thermocouples with 
stainless steel sheaths were used to monitor solution temperature 
during experiments. The reaction temperature was constantly 
monitored and controlled by a combination of water spray cool- 
ing and water bath heating. Cooling was achieved by spraying 
tap water at approximately 17°C along the neck of the flask 
when the solution temperature began to deviate from the desired 
temperature. 

Fig. 1 shows the experimental setup, excluding the water 
spray cooling. The temperature of the water bath was maintained 
constant within a variation of £0.1 °C using a thermostatic circu- 
lator. The reaction solution temperature was maintained within a 
variation of +1 °C. The NaBH; solution was also premixed with 
3.75% NaOH to suppress the hydrolysis reaction during stor- 
age (the PH of the solution was approximately 14). Two 25 mL 
flasks were used. Because of the difficulties associated with 
solid dispensing during gas-evolving experiments, the desired 
amount of ground solid catalyst in powder form was preloaded 
into the reaction flask together with a magnetic stirring bar. To 
eliminate the influence of initial wetting of the catalyst on the 
measured kinetics, 2 mL of de-ionized water was injected into 
the reaction flask to pre-wet the catalyst. The reaction flask was 


Gas burette 


Heating Flask 


Reaction|Flask ab 


Heating Mantle Water Bath 


Fig. 1. Schematic of the experimental setup. 
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preheated by submerging it for 5 min in the water bath main- 
tained at the desired reaction temperature. During this period 
the flask was stirred at 800 rpm. Then, 12 mL of 0.993% sodium 
borohydride aqueous solution with 3.75% sodium hydroxide as 
a stabilizer was injected into the heating flask. The solution was 
quickly heated from room temperature to several degrees above 
the desired reaction temperature using a small heating mantle. 
The purpose of the heating mantle is to minimize the heating time 
of the solution, which in turn minimizes the self-decomposition 
of the sodium borohydride aqueous solution during the heating 
process. In addition, using a heating mantle instead of a water 
bath enabled us to raise the solution temperature several degrees 
above the reaction temperature to account for heat losses during 
the process of transferring the solution from the heating flask to 
the reaction flask. 

Approximately 8 mL reaction solution was then transferred 
to the other pre-heated flask by a syringe to initiate the reac- 
tion. The exact amount of solution transferred was determined 
by weighing the syringe before and after injection using an ana- 
lytical balance. The maximum injected amount of solution was 
dictated by the capacity of the gas-burette, which had a mea- 
surement range of 0-250 mL. After accounting for the water 
used to wet the catalyst, the resulting solution’s initial concen- 
tration of sodium borohydride was approximately 0.8% and the 
concentration of sodium hydroxide was approximately 3.0%, 
which is commonly used for hydrolysis purposes [5,11]. A type- 
T thermocouple with a stainless steel sheath was inserted into 
the reaction flask to monitor the reaction temperature continu- 
ously. The Hz generated was measured in a gas burette by the 
volume of water it displaced. 


3. Experiments 


To establish a standard test procedure, a series of experi- 
ments was conducted to study the effects of catalyst support, 
size of catalyst particles, catalyst loading, and reactor mixing 
conditions before beginning kinetic measurements at different 
temperatures. The following section discusses the experimental 
results from these tests. 


3.1. Effects of substrate 


Three commercially available Ru catalysts from Johnson 
Matthey: 2 wt.% Ru on alumina pellets (3 mm in diameter and 
5mm in length, specific area ~200 m? g~!), 3 wt.% Ru on 
carbon extrudate (2mm in diameter and 3 mm in length, spe- 
cific area ~1000 m? g7!) and 0.5 wt.% Ru on carbon granules 
(1.7-5 mm, specific area ~1000 m? g~'), were tested for stabil- 
ity and activity at room temperature (the constant temperature 
bath was not used for these three tests). The same total mass 
of catalyst metal of each type was used in the tests. We found 
that the Ru on alumina substrate broke apart severely and crum- 
bled during the experiment. The reason could be that alumina 
is slightly acidic and reacts with the highly basic sodium boro- 
hydride solution because of the presence of sodium hydroxide 
as a stabilizer. The Ru on carbon granules also experienced 
significant break-up during experiments, perhaps due to the 


Table 1 
Sieve set and catalyst size 


Catalyst Mesh range Range of catalyst Average catalyst 
number size (um) size (um) 

1 #30-#40 425-600 512.5 

2 #50-#60 250-300 275 

3 #100-#140 106-150 128 

4 #200—#230 63-75 68 

5 #325—#400 38—45 41.5 

6 #400—#450 32-38 35 

7 #450—#500 25-32 28.5 

8 #500—#635 20-25 22.5 


irregular shapes of the particles. In comparison, 3% Ru on car- 
bon extrudate exhibited the highest integrity during experiments; 
therefore, it was used in subsequent experiments. 


3.2. Effects of catalyst particle sizes 


To study the effect of catalyst particle size, 3% Ru on 2 mm 
carbon extrudate was ground and sieved into different sizes as 
described in Ref. [12]. The twelve different sieves were stacked 
on top of each other, and the average diameters of the supported 
catalyst particles trapped in each sieve are assumed to match 
the average pore sizes of the two adjacent sieves. The sieve set 
used and the average catalyst size in each category are listed in 
Table 1. 

The effects of different catalyst particle sizes were studied at 
25°C and at 55°C. Fig. 2 shows that at 25°C, hydrogen evolu- 
tion curves for 41.5 um and 35 um overlap each other, implying 
that the reaction rate did not change with a decrease of catalyst 
particle size from 41.5 um to 35 um. Fig. 3 shows similar results 
at 55°C for catalyst particle sizes from 41.5 um to 28 um. The 
results indicate that the 35 zm catalyst was free from internal 
diffusion for temperatures from 25 °C to 55 °C, and as a result, 
35 wm catalyst particles were used to study that temperature 
range. For temperatures from 60°C to 85°C, 28.5 um catalyst 
particles were used to minimize internal diffusion effects. The 
reason that the 28.5 um catalyst was not used from 25°C to 
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Fig. 2. Hydrogen generation as a function of time for different catalyst particle 
sizes (25°C, 200 mg catalyst). 
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Fig. 3. Hydrogen generation as a function of time for different catalyst particle 
sizes (55 °C, 21.5 mg catalyst). 


55 °C was the limited amount of the catalyst available from the 
grinding process. 


3.3. Effects of stirring speed 


To obtain intrinsic kinetic data independent of external diffu- 
sion effects, four tests were conducted at room temperature using 
four different stirring speeds, 750 rpm, 1000 rpm, 1250 rpm and 
1500 rpm, and hydrogen evolution results are illustrated in Fig. 4. 
The curves for 750 rpm and 1500 rpm lay on top of each other, as 
do the curves for 1000 rpm and 1250 rpm. The small gap between 
these two pairs of curves is indicative of the experimental error. 
Therefore, the results suggest that the effect of mass diffusion 
across the liquid-phase boundary to the catalyst surface was 
eliminated when the stirring speed was at 750 rpm. For higher 
temperatures, liquid-phase diffusivity increases. As a result, a 
lower stirring speed is adequate to achieve complete mixing 
at higher temperatures. Therefore, the 750 rpm stirring speed 
is considered to be sufficient for higher temperatures as well. 
For all the tests reported hereafter, the stirring speed used was 
1500 rpm unless otherwise stated. We note that when the solution 
was not stirred, the measured reaction rate was much smaller. 
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Fig. 4. Hydrogen generation as a function of time at different stirring speeds 
(25°C, 200 mg, 41.5 um size catalyst). 


Table 2 
Catalyst amounts used at different temperatures 


Temperature (°C) Weight (mg) Catalyst size (wm) 


25.0 200.0 35 
31.1 137.6 35 
34.9 87.9 35 
40.7 52.7 35 
45.5 34.7 35 
50.9 21.8 35 
56.3 21.3 35, 28.5 
60.0 13 28.5 
64.0 10.6 28.5 
71.5 8.9 28.5 
75.8 7.1 28.5 
79.9 6.2 28.5 
84.5 5.1 28.5 


3.4. Effects of temperature 


The effects of temperature were studied from 25°C to 85°C 
at intervals of approximately 5°C. The catalyst size used from 
25°C to 55°C was 35 um while the catalyst size used from 
60°C to 85°C was 28.5 um, as explained above. Because of 
the increase in the reaction rate, the amount of catalyst used 
was reduced to a level such that the reaction process could be 
monitored accurately. The amounts of catalyst used in these tests 
are listed in Table 2. 


4. Data analysis method 


To analyze data from experiments conducted at different tem- 
peratures, we must convert raw data for hydrogen generation as 
a function of time to reaction rate versus sodium borohydride 
concentration or sodium borohydride concentration versus time. 
Assuming complete conversion of the NaBH;4 at the end of each 
experiment, the concentration at each time point based on the 
initial concentration, hydrogen generated, and the total hydrogen 
collected was calculated as follows: 


1— V 
Ca = Cao (=) (2) 
Vu, max 


Total hydrogen collected at the end of the experiment was also 
used to verify complete consumption of the borohydride. The 
small differences between measured and expected hydrogen pro- 
duction and consistency in total hydrogen generated (less than 
+5% variation in asymptotic values in Figs. 2—4) over the broad 
parameter ranges studied serve to validate this approach. Initially 
a differential method was used to analyze the data in which the 
reaction rate was calculated based on the ratio of the difference 
between two successive generated hydrogen readings to the time 
interval [12], but the resulting reaction rate versus concentration 
exhibited significant scatter because of experimental errors asso- 
ciated with each reading. Instead, we used an integral method 
[13] to analyze the data. To use the integral method for batch 
reactor data, a kinetic model (such as zero-order, first-order and 
nth-order) must be specified. The mathematical derivations for 
the integration format based on four different models are given 
in the subsequent sections. 
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4.1. Zero-order kinetics 


For a batch reactor with a volume V and a catalyst mass Meat, 
the reaction rate per unit volume based on zero-order kinetics 
can be described as: 

dCa Meatko 


= = 3 
rA dr V 6) 


Note that kọ has units of mol (masscat time)” !. Integrating Eq. 
(3) and noting that Ca cannot be negative, we obtain: 
aia- SO) orrs AE (4) 
V Meatko 

Therefore, a plot of Cao — Ca as a function of time should give 
a straight line, and the slope of the line can be used to calculate 
the zero-order rate constant ko. In principle, the line should pass 
through the origin only if the reaction is zero-order during an 
initial period. We discuss initiation by sorption of the NaBH4 
on the surface of the catalyst as a possible cause of a non-zero 
intercept in Section 5.1. 


4.2. First-order kinetics 


For a batch reactor with a volume V and a catalyst mass meat, 
the reaction rate per unit volume based on first-order kinetics 
can be described as: 

dC, Meatky Ca 


= E 5 
"A dr V ©) 


The units of kı are vol (massScat time)! . Separating and integrat- 
ing, we obtain: 


Ca dC aki ff 
| A _ Meat | dt (6) 
Cao CA V Jo 
atk 
ry ec ere (7) 
CA V 


Therefore, a plot of In(CĘo/Ca) as a function of time should give 
a straight line, and the slope of the line can be used to calculate 
the first-order rate constant ky. 


4.3. nth-order kinetics 


For a batch reactor with a volume V and a catalyst mass Meat, 
the reaction rate per unit volume based on nth-order kinetics can 
be described as: 

dCa Meatkn Ch 


a y (8) 


n—1 


The units of kn are vol” (mol” ~ ! masscat time)~!. Separating and 


integrating, we obtain: 


Ca dc atka f 
I = o Meatkn I dt (9) 
Cao CA V 0 
1 1— 1— m atk. 
rO C=- y OAD (10) 


Therefore, a plot of (1/(1 — n)(CAQ” -C a ) as a function of 
time should give a straight line through the origin, and the slope 
of the line can be used to calculate nth-order rate constant kn. 


4.4. Langmuir-Hinshelwood Model 


The reaction studied here involves a liquid phase reaction on 
a catalyst surface. We assume that the reaction consists of two 
important steps. The first step is the equilibrated adsorption of 
BH, on the surface of the catalyst as shown in the following 
equation: 


BH4 +x* = BH4 -* (11) 
The surface coverage 0a of adsorbed species is given by the 
Langmuir adsorption isotherm [14]: 


_ KCA 
~ 1+KCa 


The second step is the reaction of the adsorbed species to form 
hydrogen as: 


BH,” * +2H20 —> BO,7. x +4H2-* (13) 


Oa (12) 


Because water is abundant, the reaction is assumed to be propor- 
tional to the quantity of adsorbed NaBHy molecules 0,4. Then, 
for a reactor with a volume V and a catalyst mass meat, the 
reaction rate per unit volume can be described as: 
MeatkLOa 

=, =—_ e 14 

y a4 
Here, kų has units mol (masscat time)—!. Combining Eqs. (12) 
and (14): 


dCa Meat KCA 
ra = = 


dt V “1KO 
Eq. (15), representing the combination of Langmuir adsorption 
isotherm and kinetics that are first-order in the adsorbed species, 
is called the Langmuir-Hinshelwood kinetic model [14]. 
The adsorption coefficient K normally varies with tempera- 
ture according to: 


A Hads 
K = A exp (25s) (16) 


(15) 


That is, K decreases with increasing temperature because the 
enthalpy of adsorption is almost always negative. As a result, at 
sufficiently high temperatures, KC, is much less than unity and 
the reaction kinetics shown in Eq. (15) become first-order. At 
lower temperatures, K becomes larger. As a result KCA is much 
greater than unity, and the reaction kinetics shown in Eq. (15) 
become zero-order. To facilitate data analysis, Eq. (16) can also 
be written as: 


AH, A Hads 
K = Koexp ads | ads (17) 
RTo RT 
Separating and integrating Eq. (15), we obtain: 
CA 1+ KC t kymeg 
-f tiac = f KLMcat 4, (18) 
Cao KCA 0 V 
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(Cao 


1 Cao KLMeat 
C l = t 19 
A) + x” ( a) V (19) 


Plotting (Cao — Ca) + (1/K)ln(Cao/Ca) as a function of time 
should give a straight line, and the slope of the line (ksiope) can 
be used to calculate kg as follows: 
kLMcz 

kstope = = (20) 
We note that the left side of Eq. (19) contains the left side of 
Eq. (4) of the zero-order model, the left side of Eq. (7) of the 
first-order model, and the adsorption coefficient K. Therefore, 
this model can be considered to be a combination of zero-order 
and first-order models, with the adsorption coefficient K deter- 
mining which term is more important. The use of the adsorption 
coefficient K enables the model to capture zero-order behavior at 
low temperatures and first-order behavior at high temperatures. 
Note, however, that the value of K must be assumed in order to 
prepare the linear plot implied by Eq. (19). 

Because the volume of the batch reactor was always 10 mL, a 
normalized reaction rate constant ky” = (kL/V) = (Kslope/Mcat) 
was used to simplify data analysis. 


4.5. Application of the kinetic models 


Because an initial lag in the reaction consistently occurred, 
some data were neglected in curve-fitting the linearized models. 
At the end of most experiments, the reaction began to deviate 
from zero-order kinetics because of the limited reactant avail- 
able, and some of the data at the end of the experiments were 
also excluded. In the analysis, we first tried both zero-order and 
first-order models to fit the data, but neither of them applied well 
to the whole temperature range from 25°C to 85°C. The nth- 
order model fit data from 25 °C to 85 °C fairly well but required 
the reaction order to increase with temperature. These results 
led us to the Langmuir-Hinshelwood model as a means to fit the 
kinetic behavior over the full temperature range. 

To use Eqs. (19) and (16) in fitting the experimental data, 
two variables in Eq. (17) must be estimated first. Ko was 
determined by manually fitting the data at 25°C where T= Tọ 
with a wide range of Ko values. Through visual observation, 
we found that only values near Ko =0.22L_mmol7! matched 
the linearity of the data, as shown in Fig. 5. The uncertainty 
range for Ko is from 0.12 L mmol! to 0.50 L mmol7!, i.e. this 
range of Ko values preserves the linearity. As also shown in 
Fig. 5, values of K=0.05L mmol! and K=1Lmmol7! give 
unacceptable fits. The enthalpy of adsorption was determined 
similarly by manually fitting the 85°C data. It was found that 
AHads = —35 kJ mol! produced a good fit to the linear portion 
of the data, as shown in Fig. 5. The range of reasonable values 
for AHags is from —30kJ mol~! to —40kJ mol™!, and Fig. 5 
shows that the fits with AHads = —20 kJ mol! or —60 kJ mol! 
are clearly unacceptable. With the assumed Ko and AAgqs val- 
ues, the Langmuir model fits the data at 55°C very well, as also 
shown in Fig. 5. 

To evaluate the validity of Ko and AHaqs obtained above, 
we turn to adsorption theory. The adsorption coefficient can 


be related to the entropy change and the heat of adsorption as 
[15]: 


AS° A Hads A Hads A Hads 
— = Ko exp = 
R RT RTo RT 


- 10° A Hads 
= 0.22 exp ELL y exp | — = 
8.31 x 298 RT 


k= exp 


—15.6 x 8.31 A Hads 
= exp 
R RT 
—130.0 A Fads 
= _ 21 
exp ( - ar ) (21) 
Therefore: 
AS° = —130kJkmol7! K7! (22) 


The fact that AS° <0 indicates that the total disorder in the sys- 
tem decreases in the adsorption process, and is thus consistent 
with expectations. 


4.6. Validation of data analysis method 


To compare the reaction rates from two different tests with 
different amounts of catalyst, the slope obtained in Fig. 5 was 
normalized by the amount of catalyst used to obtain normalized 
reaction rate constant, kj’. To verify the validity of this method, 
two tests using different amounts of catalyst were conducted at 
room temperature. One used 150 mg of catalyst while the other 
used 200 mg of catalyst. The difference between the measured 
ki’ values is within 3%, further supporting the validity of the 
data analysis method used in this work. 


5. Results and discussion 
5.1. Adsorption of reactant on the catalyst surface 
At lower temperatures, the time needed for the initial H2 


bubble to appear was much longer than that at higher tempera- 
tures. In the curve fitting, these initial data points were ignored. 


© Data-25°C, Ko=0.22 Limmol 
—*—CF data-25°C 

A Data-55°C, Ko=0.22, Hads=-35 
-*- CF data-55°C 

© Data-85°C, Haas =-35 kJ/mol 
oY =172x -15 --*--CF data-85°C 


J Rĉ?= 0.998 


x Ko=0.05 L/mmol 
X Ko=1.00 L/mmol 
O Hags=-20 kJ/mol 
m Haas=-60 kJ/mol 


(Cao-Ca)+1/KIn(Cao/Ca) 


0o 2 4 6 8 10 12 14 46 48 
Time t (min) 


Fig. 5. Data analysis based on Langmuir-Hinshelwood model. CF data means 
data used in curve-fitting process, ignoring the initial data points. 
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Fig. 6. The inverse of startup time as a function of the inverse of temperature. 


We assume that this initial period occurred because the catalyst 
surface had not yet reached the adsorption equilibrium that is 
assumed in the Langmuir-Hinshelwood model discussed above. 
We can calculate the induction or sorption time by extrapolating 
the linear fit to its intercept with the time axis. The inverse of the 
sorption time is related to the sorption rate. The resulting esti- 
mates of the sorption rates are plotted in Fig. 6. It is interesting 
to note that this measure of the sorption rate varies exponen- 
tially with inverse temperature, which is expected Arrhenius 
behavior. 


5.2. Effects of stirring speed 


The experimental data in Fig. 4 were analyzed using the 
Langmuir-Hinshelwood model, and results are compiled in 
Table 3. The normalized reaction rates match the theory 
well, within +3%, indicating that external diffusion effects 
were insignificant when stirring speed varied from 750 rpm to 
1500 rpm. This result also validates the conclusions in Section 
3.3; 


5.3. Effects of catalyst size 


The Langmuir-Hinshelwood model was used to analyze data 
in Figs. 2 and 3, and the results are compiled in Tables 4 and 5. 
At 25°C, the normalized reaction rate constant kj’ in general 
increased with decreasing catalyst size, doubling as the catalyst 
size decreased from 512.5 um to 35 wm. At 55 °C, the normal- 
ized reaction rate constant kj’ also increased with decreasing 
catalyst size as expected because the true rate constant generally 
increases faster with temperature than the effective diffusiv- 
ity, thereby causing the effectiveness factor to decrease at a 
given particle size and the effect of particle size to be more 
marked at higher temperatures. The effective value of ky’ 
increased almost 20-fold when the catalyst size decreased from 
512.5 um to 35 um for reactions at 55°C. The small varia- 
tion in effective ky” for catalyst diameters of 28.5-41.5 um 
suggests that the rate is independent of particle size in this 


range. 


5.4. Effects of temperature and activation energy 


Experimental data at different temperatures ranging from 
25 °C to 85 °C were analyzed using the Langmuir-Hinshelwood 
model, and the results appear in Table 6. All correlation 
coefficients were larger than 0.99, implying the Langmuir- 
Hinshelwood model captures the reaction process at different 
temperatures very well. We note that normalized reaction rate 
constant, ky’ , increased more than 130-fold when the tempera- 
ture increased from 25 °C to 85°C. 

The reaction rate constant, ki’ , at different temperatures can 
be plotted on an Arrhenius plot to calculate the activation energy 
of the catalyzed reaction as shown in Fig. 7. A good fit occurs 
with a correlation coefficient of 0.975, supporting the assump- 
tion that the 28.5 um catalyst size is sufficient to avoid significant 
diffusion influences in the 55—85 °C temperature range. We note, 
however, that a careful examination of Fig. 7 reveals that the ini- 
tial slope of the curve is higher than that of the later part of the 
curve, perhaps due to the onset of diffusion effects at higher tem- 
peratures or to experimental error. The high-temperature data 
exhibit significant scatter because of errors associated with tem- 
perature control and measurement of the small weight of catalyst 
used to accommodate the higher activity. Taking data only from 
25°C to 55°C, we calculate the upper bound of the activation 
energy to be 83.5 kJ mol~!, which is about 20% higher than the 
calculated value of 66.9 kJ mol! based on data from 25°C to 
85°C. 

From Fig. 7, the relation between the normalized reaction 
rate and temperature can be written as: 


8054.6 
Tk 


ki” (mmol L~! min! mg cat.~!) = exp (- + 25.601) 


(23) 


The above equation can be used to calculate the rate constant k, 
and to convert it to standard units as: 
8054.6 
) (24) 
T, 


K 


kı (kmols~! kg cat.~!) = exp(16.903) exp (- 


Eqs. (15), (17) and (24) represent the final form of the present 
kinetic model developed for reactor design and modeling. 
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Fig. 7. Arrhenius plot for sodium borohydride hydrolysis on ruthenium catalyst. 
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Table 3 

Effects of stirring speed on kj” (25°C) 

Stirring speed (rpm) Weight (mg) Rate equations? [y (mmol L~!), x (min)] Correlation coefficient, R? ki” (mmol L! min“! mg!) 
750 200.1 y=30.3x— 155 0.997 0.151 

1000 198.6 y=30.0x— 175 0.998 0.151 

1250 200.1 y=29.7x— 170 0.999 0.148 

1500 199.4 y=30.9x — 168 0.998 0.155 


à Equation y=ax + b comes from linear fitting of the experimental data based on Eq. (19). K in Eq. (19) was calculated based on Eq. (17). The non-zero intercept 
b is due to the initiation process explained in Section 5.1. The slope a is equal to kmcat/V based on Eq. (20). 


Table 4 
Effects of different catalyst particle sizes on kj’ (25°C) 
Catalyst size (wm) Weight (mg) Rate equations* (mmol L~! vs. min) Correlation coefficient, R? ki” (mmol L! min“! mg!) 
512.5 200.3 y=14.3x — 106 0.999 0.071 
275 199.3 y= 12.6x — 63.8 0.999 0.063 
128 199.5 y=21.5x— 110 0.999 0.108 
68 199.3 y=26.9x — 140 0.997 0.135 
68 199.5 y=27.9x — 147 0.999 0.140 
41.5 199.4 y=32.5x— 178 0.998 0.162 
35 199.6 y=31.6x— 179 0.997 0.158 


à Refer to the note in Table 3. 


Table 5 
Effects of different catalyst particle sizes on kj’ (55°C) 
Catalyst size (wm) Weight (mg) Rate equations* (mmol L~! vs. min) Correlation coefficient, R? ky" (mmol L! min“! mg!) 
512.5 21.2 y=4.4x+46 0.906 0.206 
275 21.5 y=15.5x+5 0.989 0.723 
128 21 y=30.4x — 16 0.993 1.447 
68 21.3 y=54.7x — 48 0.999 2.568 
41.5 21.5 y=79.3x — 68 0.999 3.688 
35 21.3 y=80.1x— 68 0.999 3.759 
28.5 21.3 y=86.9x — 72 0.997 4.078 
è Refer to the note in Table 3. 
Table 6 
Effects of temperature 
Temperature (°C) Weight* (mg) Rate equations? (mmol L~! vs. min) Correlation coefficient, R? ki” (mmol L! min“! mg!) 
25.0 200 y=32.5x— 178 0.998 0.162 
31.1 137.6 y=46.5x — 184 0.999 0.338 
34.9 87.9 y=51.4x— 165 0.999 0.585 
40.7 52.7 y=60.4x — 160 0.998 1.147 
45.5 34.7 y=60.5x — 134 0.999 1.745 
50.9 21.8 y=53.6x — 80 0.999 2.461 
56.3 21.3 y=86.9x — 72 0.997 4.078 
60.0 13 y=66.9x — 41 0.998 5.146 
64.0 10.6 y=63.2x —21 0.998 5.966 
71.5 8.9 y=81.5x— 12 0.997 9.157 
75.8 7.1 y=60.1x+15 0.994 8.461 
79.9 6.2 y=78.0x+4 0.993 12.573 
84.5 5.1 y=109x —17 0.998 21.274 


è Refer to Table 2 on the information of catalyst size. 
> Refer to the note in Table 3. 
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Table 7 
Effects of different catalyst loading’ 


Catalyst loading Weight (mg) Rate equations? (mmol L~! vs. min) Correlation coefficient, R? ki" (mmol L! min“! mg!) 
3% 199.4 y=32.5x— 178 0.998 0.163 
4.41% 136.0 y=67.2x — 223 0.995 0.336 
10% 59.8 y=107x — 254 0.997 0.539 
5% wet 238.9 y=157x— 168 0.996 0.789 


è The sizes of all catalysts were 41.5 um except for 5% wet, whose size was not available because the catalyst all clusters together due to water. 
b ki" was calculated based on the equivalent weight of 3% Ru on carbon to facilitate comparison between different catalyst. 


€ Refer to the note in Table 3. 


275+4------ —+- Size 41.5um ground 3% Ru on carbon extrudate |-------- 
2501------ = Size 41.5um 4.41% Ru on carbon powder 
-® Size 41.5um 10% Ru on carbon powder 

2254------ ~@ 5% Ru on carbon powder, 50% water wet 
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Fig. 8. Hydrogen generation as a function of time for different catalyst loadings 
(25 °C, 6.0 mg net Ru). 


5.5. Effects of different catalyst loadings 


To study the influence of the catalyst loading, four tests were 
conducted at room temperature. The first test used ground 3% Ru 
on carbon extrudate. The second test used 4.41% Ru on activated 
carbon powder, while the third used 10% Ru on activated carbon 
powder. The last test used 5% Ru on activated carbon powder, 
nominally 50% water wet (i.e. the final catalyst contains 50% 
of water by weight). The net amount of Ru in each test was 
maintained constant, and the size of catalyst used was 41.5 wm 
except for the 5% catalyst, whose size was not available because 
the catalyst particles clustered together. Fig. 8 shows hydrogen 
generation as a function of time for different catalyst loadings, 
and Table 7 contains the results of the analysis. In principle, 
the differences in observed rates could derive from the highly 
exothermic nature of the reaction, but the spread of k{” by a 
factor of 4.8 would require a temperature rise of approximately 
18°C, which is unlikely considering the small catalyst particle 
size and the liquid medium available for heat transfer. The vari- 
ations in rate imply that the weight of Ru alone is not sufficient 
to characterize the catalysts. Ru particle size, Ru distribution on 
the particle surface as opposed to within the volume, and the 
influence of the carbon surface hydrophobicity on accessibil- 
ity of the pore structure to the reacting liquid are some of the 
properties could all differentiate the catalysts. The differences 
in catalytic efficiency of the Ru, shown in Fig. 8 and Table 7, 
indicate that further optimization of the process may be possible 


and that more detailed catalyst characterization is warranted in 
future work. 


6. Conclusions 


A detailed kinetic study of sodium borohydride hydrolysis 
on Ru-on-carbon catalyst has been conducted. Experimental 
data show zero-order kinetic behavior at low temperatures and 
first-order kinetic behavior at high temperatures. Based on exper- 
imental observations, a unified kinetic model was developed. 
The results indicate that the reaction consists of two impor- 
tant steps: the adsorption of sodium borohydride on the surface 
of the catalyst and the reaction of the adsorbed species on 
the catalyst surface. These steps are combined in the classical 
Langmuir-Hinshelwood model that shows that the adsorption 
step is responsible for variation of kinetic order of the reac- 
tion from low to high temperatures. Furthermore, reaction rate 
constants at different temperatures were determined from exper- 
imental data, and the activation energy of the catalyzed process 
was determined from an Arrhenius plot. 
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